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Abstract 

Using the so(2, 1) Lie algebra and the Baker, Campbell and Hausdorff 
formulas, the Green's function for the class of the confluent Natanzon po- 
tentials is constructed straightforwardly. The bound state energy spectrum 
is then determined. Eventually, the three-dimensional harmonic potential, 
the three-dimensional Coulomb potential and the Morse potential may all be 
considered as particular cases. 
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In this paper, we shall want to present an algebraic treatment of the bound 
state problem for the class of the confluent Natanzon potentials defined by 



V{r) = h 



R i?2 4 ^3 ' 

where A = a\—Aa2Co and gi, g2, ci, o"2, cq and r] are dimensionless parameters. 
The function h = h[r) is defined implicitly by the differential equation 

dh ■ 2h 

dr .[^ 



where 

R{r) = a2h'^ + (Tih + Co. (3) 

These potentials, originally introduced by Natanzon, are interesting in the 
sense that the Schrodinger equation can be reduced to the confluent hyperge- 
ometric form and they include, as special cases, the radial harmonic oscillator 
for (72 = Co = 0, the radial Coulomb potential for ai = Cq = and the Morse 
potential for (T2 = cti = 0. These confluent Natanzon potentials have been 
recently the object of several studies. More specifically, NatanzonB has ob- 
tained, in a Schrodinger equation approach, the discrete energy spectrum, the 
unnormalized wave functions and has derived the Green's function. Cooper 
and all have used the supersymmetric quantum mechanics as an algebraic 
method to discuss these potentials. Furthermore, the spectrum of these po- 
tentials has been re-derived by Cordero and Salamoi , still in a Schrodinger 
approach, by making use a particular realization of the so{2, 1) algebraic 
method which is called spectrum generating algebra (SGA)i . ' Finally, this 
class of the confluent Natanzon potentials has been solved very recently by 
Groschei in the framework of a Feynman path integral approach. 

The Milshtein and Strakhovenko variant (MS) of the so(2, 1) algebraic 
approach! provides an elegant and straightforward way of solving the problem 
for V{r). As we shall see, this approach consists of constructing the Green's 
function by expressing the resolvent operator for this system in the form of a 
linear combination of an appropriate realization of 5*0(2, 1) generators. This 
form is then written in the Schwinger's integral representation and the action 
of the resolvent operator is obtained by means of two Baker, Campbell and 
Hausdorff formulas. 

Let G{r, r'; E) be the Green's function associated to the potentials U{r) = 
^V{r), which is solution of the differential equation 

(H - E)G(r,r';E) = -6(r -r'), (4) 
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where 

is the Hamihonian of the system and E its energy. 

In terms of the function h{r) and by taking into account the relation (|]), 
we can write the resolvent operator H — E as follows: 

2m \ dr^ ^ \ h h 

^ \g2 - (J2e)h^ + {gi - cTie)h + {t] - 1 - coe) } , (6) 



with E = ^e. 

2m 



The point canonical transformation ^ = h{r) leads to 



H-E)Gi^,<;,E)^\m-a (7) 



where 



G{i,i'-E) = {r)h{r')G{ry-E). 



Then, the dynamics of the physical system is governed by the new Hamilto- 
nian 



/_ 7] -I- Coe , g2 - o-ae ^ 

2m 1^ rff 4e 



H=^A-^^, + "-\;^'^ +'-^^]. (9) 



and 

E = -^{g,-a,e) (10) 

is the pseudo-energy. 

To see that one can introduce the three following generators 

(11) 

satisfying the so(2, 1) Lie algebraB 

[Ti, T2] = -«Ti, [T2, Ts] = -zTa, [Ti, T3] = -zTa, (12) 
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it is convenient to perform the change of variable = x'^ and redefine the 
Green's function as 



G{^, E) = V^G{x, x'; E). 
Under this change the differential equation (|^) becomes 



n 



(Ti(x) + 2h^u^Ts{x) - e) G{x,x- E) = -6{ 



X — X 



where 



(13) 

(14) 
(15) 



with M = Am. 

Being a linear combination of the generators Tj, the operator H shows, as 
expected, a dynamical symmetry so{2, 1). By using the Schwinger's integral 
representation! , the solution of the differential equation (|l^) can be written 
as follows: 



where 



G(x,x';E)= / dS exp 
Jo 



iS 



[E + i0)s K{x,x';S), 



(16) 



K{x, x'; S) = exp |-y [Ti(x) + 2h^u^T3{x)\ j 6{x - x'). (17) 

The calculation of this kernel is based upon the use of two Baker, Campbell- 
and Hausdorff formulas^ 



iS 

exp i -- 



Ti + 2%^uj'^T'3 \ = exp(— iaTa) exp(— ^6X2) exp(— icTi), (18) 



where 



(19) 



a = 2hujtaiiiujS), 6 = 2 In [cos(c<jS')] , c = — tan(tt;S'), 



and 



exp(— iaTs) exp(— i/3T2) exp(— z7Ti) = exp(— icTi) exp(rT3), 
with 

i'T ^ iTC\ c 

0^ = - /? = 21n 1-— , 7 



1 - iTc/2'- 



1 - iTc/2' 



(20) 



(21) 



It is to be noted that formulas (|T8|) and ( pOf ) can be easily checked using the 
finite representation of so{2, 1) Lie algebra where the Tj operators are defined 
in terms of a, Pauli matrices 



Ti 



rr - 

i9 — —■ 



Ox + 102 



Here, we also have to use the Laplace transform of the Dirac distribution 



(22) 



b{x — x') 



2h^ 2m 



ioo+S 



drexp 



M 
4^ 



(x^ - x'^)r 



6 ^ 0, 



(23) 



in order to obtain a manageable result as follows: 



exp(— i7Ti)a;'' 



1 - i'jTi 



2! 



-i'jTi 



x^" = x^. 



(24) 



Using relations (pSj), (|T8|) and (|20|), the kernel ( p!?] ) can now be written 

K(x,x';S) = — oO;' exp(— zaTs) exp(— z6T9) / dr 

2n 2zn J-ioo+s 

f M ,0 \ , ^ , ( ^ 'l\ u 

X exp I lyX T exp(— 2cTi) exp — ^x r x^ 

\ Ah J \Ah / 

= — 2^ ''exp(-mr3) exp(-i6T2)— / dr 

2h 2t7T J-ioo+5 



X exp 



M 



Ah' 



X T ) exp(— icTi) exp(rT3)x'' 



M ,1 1 fioo+S 

—^x expi^-iaTs) exp{-ibT2)— / dr 
2h 2m J-ioo+5 



X exp 
2/1^ 



M 



x' T I exp(— iaTa) exp(— 'i/9T2) exp(— i7Ti)x^ 



■ exp(— mTa) exp{—ibT2)x^ 



2m 



ioo+5 

X / dr 

joo+<5 



M+7 



2h^ 



tX 



' exp(— mTa) exp{—ibT2)x^ exp 



exp (-^rj exp (^7^ 



X / (ir 

- joo+(5 



where formulaB0 



(1-^) 



(25a) 



exp(— i/9T2)/(x) = exp 



4)/(e-|.) 



(26) 



has also been used. 

The integral can be calculated thanks to the residue theorem after the 
^^P (|^ T--|-2i/c ) series has been effected. Hence we obtain 



K{x,x';S) 



M 

— 2~ 6xp ( —iaT^ ) exp ( —ibT2 ) v xx' exp 
ih c 



iM 



2rc 



_ 2 , 
2 y i-c I 



2 



Mxx' 



ih c 



M 



exp 



iaM 
7rX~ I exp 



X exp 
Muj 



iM 



2h^c 



[e-'x^ + x") 
iMu) 



ih sin(tt;S') 



' XX exp 



2h 



'xx' 
Mxx' 



M-2 \ •f;2 i 

\th ce^ . 



(x^ + x'^) cot(u;S') 



2 



Mujxx' 
ih sm(u) S) ^ 



where 



(27) 



(2^ 



and /,,_ 1 (x) is the modified Bessel function. 

After inserting ( P7D into (JT^) and performing the integration over the 
time variable 5* with the help of the following formula (see Gradshtein and 
Ryzhik0 , p. 729 , eq. (6.669.4)) 



Jo SI 



-2pq 



sinh q 



exp 



— -(x + y) cothg 



'27 



xy 



sinh q 



r(27 + l)yxy 



M_p,^(x)iy_p,^(y), 



(29) 



valid for i?e (^p + 7 + ^ > 0, -Re(7) > and y > x, where M_p^^(x) and 
iy_p,^(y) are the Whittaker functions, the Green's function ([T6|) becomes 



(^(x, x'; E) 



r(p+f + i) 



2i(i;r(/i + i; 



'XX 



-f^' 2 4 



x'^)W 



.Muj 



'''2 4 



-X 



(30) 



where p 



■-f- and X > x'. 
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Eventually, by taking into account eqs. (pOD, (0), (||), (0) and the suc- 
cessive transformations C, = h{r) = x^, the Green's function associated with 
the confluent Natanzon potentials may be written: 



G(r, r';E) 



r - 



2hu 



ii _|_ i 
2 4 



,Muj 



R{r)R{r') 
h{r)h{r') 



^'2 4 



IT 



/i(r')) 



^'2 4 



/^(r)), 



(31) 



and the bound state energy eigenvalues = can be obtained from the 

p°i«s + f + in®, 



(32) 



E li I 



Upon substitution of (pUj), (|TB|) and ( ^5] ) into (|52|), we see then that 
91 - oxt 



2V5'2 - cr2e 



+ - coe = -(2n + l) 



(33) 



IS an equation of fourth degree m e„ which can be solved after some simple 
manipulationsHj . 

Many particular cases of the potentials (|I|) arise when some of its param- 
eters vanish. In the following table, we give the interesting ones, namely, 
the radial harmonic oscillator, the radial Kepler- Coulomb potential and the 
Morse potential 



0-1 


0-2 


Co 


hir) 


V{r) 












91 _|_ 92 ^2 _|_ ^-J 
ai 'a^ ~r J.2 












92 1 gi 1 

o"2 2y^r' ' 4r2 












92 g~2r/co _|_ 91 g-r/v^cS" i V, 
cn cn cn 



In conclusion, we have shown that the class of confluent Natanzon po- 
tentials can be described by a particular realization of the so(2, 1) which 
enables us to work out the Green's function. The realization of the so{2, 1) 
algebra used here differs substantially from that of Cordero and Salamoi . 
The discrete spectrum is identical with the results obtained by using either 
the (SGA) method! or the Feynman path integral approach! . 
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